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We show that the abundance and redshift distribution (dN/dz) of galaxy clusters in future high- 
yield cluster surveys, combined with the spatial power spectrum (P c (k)) of the same clusters, can 
place significant constraints on the evolution of the dark energy equation of state, w — w(a). 
We evaluate the expected errors on w a — —dw/da and other cosmological parameters using a 
Fisher matrix approach, and simultaneously including cluster structure evolution parameters in our 
analysis. We study three different types of forthcoming surveys that will identify clusters based 
on their X-ray emission (such as DUO, the Dark Universe Observatory), their Sunyaev-Zel'dovich 
(SZ) decrement (such as SPT, the South Pole Telescope), or their weak lensing (WL) shear (such 
as LSST, the Large Synoptic Survey Telescope). We find that combining the cluster abundance 
and power spectrum significantly enhances constraints from either method alone. We show that the 
weak- lensing survey can deliver a constraint as tight as Am„ ~ 0.1 on the evolution of the dark 
energy equation of state, and that the X-ray and SZ surveys each yield Aw a ~ 0.4 separately, or 
Aw a ~ 0.2 when these two surveys are combined. For the X-ray and SZ surveys, constraints on 
dark energy parameters are improved by a factor of two by combining the cluster data with cosmic 
microwave background (CMB) anisotropy measurements by Planck, but degrade by a factor of two 
if the survey is required to solve simultaneously for cosmological and cluster structure evolution 
parameters. The constraint on w a from the weak lensing survey is improved by ~ 25% with the 
addition of Planck data. 



I. INTRODUCTION 

It has long been realized that clusters of galaxies pro- 
vide a useful probe of fundamental cosmological param- 
eters. The formation of the massive dark matter po- 
tential wells is dictated by simple gravitational physics, 
and the abundance and redshift distribution of clusters 
(dN/dz) should be determined by the geometry of the 
universe and the power spectrum of initial density fluc- 
tuations. Early studies of nearby clusters used this re- 
lation to constrain the amplitude eg of thepower spec- 
trum (e.g., [EQl)- Subsequent works (e.g., fcjElOl) have 
shown that the redshift-cvolution of the observed cluster 
abundance places useful constrains on the matter density 
parameter J7 m . 

The next generation of surveys, utilizing the Sunyaev- 
Zel'dovich effect (SZE), X-ray flux, or weak lensing sig- 
natures to identify galaxy clusters, will be able to deliver 
large catalogs, containing many thousands of clusters, 
with complementary selection criteria. Such forthcoming 
datasets have rekindled a strong theoretical and exper- 
imental interest in galaxy clusters. In particular, Wang 
& Steinhardt p first argued that the cluster abundance 
can be used to probe the properties of dark energy, and 
Haiman, Mohr & Holder showed that a survey with 
several thousand clusters can yield precise statistical con- 
straints on both its density (Ode) an( i its equation of 
state (w = P/p). Several subsequent recent works have 
focused on various aspects of extracting cosmological pa- 
rameters from high-yield, future surveys, such as the sta- 



tistical constraints available on curvature Qk @; assess- 
ing the impact of sample variance |j and other uncertain- 
ties [l(| on parameter estimates; and controlling such un- 
certainties by utilizing information from the shape of the 
cluster mass function dN/dM [ll|. Closest to the sub- 
ject of the present paper, Weller et al. [I^] and Weller & 
Battye [13j considered constraints on the time evolution 
of the dark energy equation of state in forthcoming SZE 
cluster surveys. 

Recent studies have elucidated the additional cosmo- 
logical information available from the spatial distribu- 
tion of galaxy clusters through a measurement of their 
three-dimensional power spectrum P c (k) [Tj| (see [l5j 
for a more general treatment of extracting cosmologi- 
cal information from redshift surveys). The power spec- 
trum contains cosmological information from the intrinsic 
shape of the transfer function |16| and from baryon fea- 
tures [13, 0, Il3 ■ The existing sample of ~ 400 nearby 
clusters in the REFLEX survey has already been used to 
derive their power spectrum; combined with the number 
counts, this has yielded constraints on erg and fi™ that 
are consistent with other recent determinations [2fj|. 

Most importantly, the cluster abundance and power 
spectrum can provide two independent powerful probes 
of cosmological parameters from a single dataset. The 
dependence of dN/dz and P c (k) on the cosmological pa- 
rameters are different. Combining the two pieces of infor- 
mation can therefore break degeneracies present in either 
method alone, and yield tighter statistical constraints. 
Furthermore, this can be used to significantly reduce 
systematic errors arising from the mass-observable rela- 



tion, making cluster surveys "self-calibrating" 21]. This 
self-calibration is especially strong when the abundance 
and power spectrum information is combined with even 
a modest follow-up mass calibration program [22] • 

In this paper, we return to the question of constraining 
the time-evolution of the dark energy equation of state. 
Specifically, we ask the question: can improved con- 
straints be obtained on the time-evolution of w — w(z) 
when the cluster counts and power spectrum are com- 
bined? We quantify the statistical constraints expected 
to be available from future samples of <; 10, 000 galaxy 
clusters. 

We study constraints from three different types of 
forthcoming cluster surveys. The proposed DUO (Dark 
Universe Observatory) [23J ] X-ray survey will be per- 
formed by an Earth-orbiting satellite consisting of seven 
telescopes that take a wide survey of the sky in soft X- 
ray bands. The SPT-like (South Pole Telescope) 0] 
SZE survey will be performed by an 8-meter precision 
submillimeter-wave telescope detecting distant galaxy 
clusters by their Sunyaev-Zel'dovich decrement. The 
LSST-like (Large Synoptic Survey Telescope) [25( survey 
will be performed by a ground-based telescope detecting 
clusters by their weak lensing shear signature on back- 
ground galaxies. 

The most important differences between the present 
paper and earlier works that have addressed the time- 
evolution of the equation of state |12l fl3 ] are that here 
(i) we simultaneously include the abundance and power 
spectrum in our analysis; (ii) in addition to the cosmo- 
logical parameters, we simultaneously include several pa- 
rameters that describe cluster structure and evolution; 
and (iii) we study three different types of forthcoming 
surveys. We also chose a different fiducial model for 
our analysis (one close to the standard "concordance" 
cosmology). Our calculational method, based on Fisher 
matrices, is, on the other hand, only a simple approxi- 
mation to the Monte-Carlo likelihood analysis performed 
in EH!. 

This paper is organized as follows. In <jnj we describe 
our basic calculational methods. In £|IIII we present our 
results for different future cluster surveys. In HTV\ we 
critically discuss our results, including their uncertain- 
ties, and summarize the implications of this work. 



II. CALCULATIONAL METHOD 

A. Simulating Cluster Surveys 

We follow the standard approach, and identify galaxy 
clusters with dark matter halos. The differential comov- 
ing number density of clusters is given by Jenkins et 
alfM 



(In 
dM 



(z,M) = QMl P^ d\na- l (M,z) 

K ' M dM 

x exp[-| lna- 1 (M, z) + 0.64| 3 ' 82 ] , (1) 



where a 2 (M, z) is the variance of the linear density field 
at redshift z, smoothed with a spherical top-hat fil- 
ter which would enclose mass M at the mean present- 
day matter density p m 56]. The Jenkins et al. mass 
function was derived from numerical simulations, and 
its self-similar form is demonstrated to be accurate to 
within ~ 15% in three widely separated cosmologies (al- 
though see |27| who find a more significant cosmology- 
dependence of the mass function) . Jenkins et al. identify 
simulated clusters using Miso, the mass enclosed within 
a spherical overdensity of 180 with respect to the mean 
matter density. However, it is customary to define the 
relation between X-ray or SZE flux and halo mass M20O) 
defined as the cluster mass enclosed within a sphere with 
mean interior overdensity of 200 relative to the critical 
density. To combine this relation with the mass function 
in Eq. QJ, we convert M200 to Miso assuming that the 
halo density profile is described by the NFW model with 
a concentration parameter of c n f w = 5 [28| . 

The spatial distribution of clusters is assumed to follow 
the spatial distribution of the dark matter halos and is 
specified by the cluster power spectrum P c (k). We follow 
Hu & Haiman [lj] and obtain P c (k) from the underlying 
mass power spectrum, P(k), modified by redshift-space 
distortions 29] 



P c (k±,h 



■^i 



A' 



II 



b 2 P(k), (2) 



(3) 



where k± and ku are the wavenumbers of the sinusoidal 
fluctuation modes transverse and parallel to the line of 
sight, respectively. The redshift-distortion parameter 8 
is defined by J23 



_ 1 d In Dg r o W 
b din a 



(4) 



where D gmw is the linear growth rate, and a is the ex- 
pansion factor normalized to unity today. The parameter 
b in Eqs. 10) and Q represents the linear bias averaged 
over all halos at redshift z: 



b(z) 



M min ( z ) 



dn(M, z) 
dM 



b(M) dM 



M min ( z ) 



dn 
dM 



-1 -1 



dM 



(5) 

where M m [ n (z) is the mass of the smallest detectable clus- 
ter, which depends on the type of survey as discussed in 
JiTlD. The bias parameter of halos of a fixed mass M, 
b(M), is assumed to be scale independent and given by 



b(M) = 1 + 



aS 2 /a 2 - 1 



2p/6 c 



1 + (aS 2 /a 2 )P 



(6) 



with a = 0.75 and p = 0.3 providing the best fits to the 
clustering measured in cosmological simulations [3(j ■ Fi- 
nally, S c represents the threshold linear overdensity corre- 
sponding to spherical collapse, whose value is 8 C = 1.686 



in an Einstein-de Sitter universe. We keep it fixed 
throughout the calculation, as it was shown to be only 
weakly dependent on cosmology and redshift in other 
models 61. 



B. Fisher Matrix Technique 



TABLE I: Parameters for the Planck Survey. 



Frequency (GHz) 


too 


143 


217 


9 C (arcmin) 


10.7 


8.0 


5.5 


<7 Ci T (nK) 


5.4 


6.0 


13.1 


(T C ,E (^K) 


— 


11.4 


26.7 


£c 


757 


1012 


1472 



The Fisher matrix formalism allows a forecasting of 
the ability of a given survey to constrain cosmological 
parameters [3l| . It gives a lower bound to the statistical 
uncertainty of each model parameter that is to be fit by 
future data. The well-known advantages of the Fisher 
matrix technique are that (i) it allows a quick estimate 
of errors in a multi-dimensional parameter space, since 
the likelihood functions do not have to be evaluated at 
each point of the multi-dimensional grid, and (ii) con- 
straints from independent data sets or methods can be 
easily combined by simply summing the individual Fisher 
matrices. 

The Fisher matrix is defined as 



F, 



d 2 C 
dpidpj 



(7) 



where C = — In L is the log-likelihood function, and 
where p^a are the various model parameters which, in 
our case, include both cosmological parameters and pa- 
rameters describing cluster structure and evolution. The 
inverse (F )y gives the best attainable covariance ma- 
trix, regardless of any specific method used to estimate 
the parameters from the data |3l| . In particular, the 
best statistical uncertainty attainable on any individual 
parameter p^, after marginalization over all other param- 
eters, is {F~ l )][ 2 . 

We construct the Fisher matrix for the redshift dis- 
tribution of the number density of galaxy clusters as |2| 



recounts 



E 



dNjdNj 1 

dp^ dp u Ni 



(8) 



where 



d 2 V 

Ni = AnAz-——(z l ) 
dzdil 



M mia (z i ) 



dn(M, Zi) 
dM 



dM 



(9) 



is the number of clusters above the detection threshold 
in each redshift bin i centered at Zj. In Eq. J5J, AJ7 
is the solid angle covered by a survey, cPV/dzdil is the 
comoving volume element, and dn/dM is the cluster mass 
function (see Eq. QJ). We sum over redshift bins of size 
Az = 0.05, between z m - m = and z max = 2.0 for the 
DUO X-ray survey and the SPT-like SZE survey, and 
between z m ; n = 0.1 and z max = 1.4 for the LSST-like 
survey, although accurate redshifts are not required for 
the dN/dz test. 



We construct the Fisher matrix for the redshift-space 
power spectrum as |14| 



ppowcr 



E 



01n(fc3fc||P c )y 01n(*i*||P c )« (VkV^ij 



dp^ dp v 2 

where P c (k) is the cluster power spectrum. The two- 
dimensional fc-space cells and the set of redshift bins 
are labeled by i and j, respectively. The factor of 
(VfeT4ff/2) _1 estimates the uncertainty (AP C /P C ) 2 in the 
measured power spectrum, including the effects of shot 
noise and cosmic variance [22J. Here V e g is the effective 
volume probed by the survey 



V eS (k) 



dV s 



n( Zj )P c (k) 



l + n{z 3 )P c (k) 



(11) 



where n is the expected average number density, and Vk 
is the cylindrical volume factor in fc-space: 



V k 



2 7 rA(fc 2 JAfc|| 
(2Vp 



(12) 



We sum over 29 2 linearly spaced fc cells from k±.\\ = 
0.005 to 0.15 Mpc , thus defining a cylinder in three- 
dimensional fc-space. We sum over redshift bins of size 
Az = 0.2, between Zmin = and z max = 2.0 for the 
DUO-like survey and the SPT-like survey, and between 
£min = 0.1 and Zmax = 1.4 for the LSST-like survey. The 
choice of this relatively wide redshift bin size is dictated 
by the need to have a sufficient number of clusters in each 
bin for an accurate determination of the power spectrum 
(N > 1, 000), as well as a wide enough Az that includes 



radial modes with fcjii w 0.005 Mpc . 

Finally, in addition to the constraints from clusters 
considered here, we construct the Fisher matrix that can 
be used to forecast cosmological parameter errors from 
the temperature and polarization anisotropy of the cos- 
mic microwave background (CMB). We have in mind a 
near-future survey such as Planck |3.1 | that will measure 
temperature and E-mode polarization auto-correlation 
(respectively, TT and EE), as well as temperature- 
polarization cross-correlation (TE). We neglect B-mode 
polarization. The full CMB Fisher matrix is then given 
by Hlli 



TTicmb 



EE 

£ X.Y 



dC 



x.. 



dpp. 



Gov" 1 iC x ,t,C Y j) 



dC Y 
dpv 



(13) 



where X,Y run over TT, EE and TE correlations. The 
covariance matrix, Cov(Cx,e, Cy,t), has elements 



Cov(Ctt,i, Ctt,i) 



-2\2 



— — (C T T,e + B Tl ) 

J- J./skv 



(2£+l)/ sky 

2 
Cov(Cee,i,Cee.i) — ln „ ; it j — {Ceej + B E 

1 J./skv 



-2\2 



(2£+l)/ sky 



Cov(Cte,£, Cte,i 



Cov(Cee,£, Cte,i) 



1 



"[Ct£,£ + 



(2£+l)/ sky ' 

(Ctt.i + B t( ){Cee,i + B E ti 



2 
Cov(CTT,e,CTE,i) = /0/; , 7Ti — Cte,i{Ctt,i + B T ,) 

{It + ij/sky 

2 2 

COV (CtT, I, CeE,1) = tn0 , 7w — C TE li (14) 



(2£+l)/ sky 



where / sky is the fraction of the sky covered. The B^/s 
and Be,iS account for experimental noise for temper- 
ature and polarization measurements, respectively, and 
are given by |^3 



b! = J2MY 



2 e -l(l+l)/l 2 c 



(15) 



matter (CDM) cosmological model, dominated by a cos- 
mological constant (A). The sensitivity of our results to 
the choice of the fiducial parameters is discussed in § I1VI 
below. The parameters are adopted from recent mea- 
surements by the Wilkinson Microwave Anisotropy Probe 
(WMAP), as summarized in Table 1 of [39j: baryon den- 
sity ttjjh 2 — 0.024, matter density D m h 2 = 0.14, dark 
energy density in units of the critical density D^e = 0.73 
(or Hubble constant Hq = lOOh km s _1 Mpc^ 1 with 
h = 0.72), with present-day normalization erg = 0.9 and 
scalar power-law slope n s = 1 of the primordial power 
spectrum. Following Linder |19J| . we parameterize the 
evolving dark energy equation of state as 



w(z) — w + w a (l - a) = w + w a 



1 + z' 



(16) 



with values in our fiducial model chosen to be wq = — 1 
and w a = 0. An alternative parameterization sometimes 
used in the literature is w(z) = w a + w z z. The errors we 
obtain here on w a should be divided by approximately 
a factor of two to obtain the corresponding errors on 
w z . This follows from Taylor-expanding Eq. (|lfj|) about 
z = 1/2, which is approximately where the sensitivity of 
cluster surveys peak. 



where the sum is over the different frequency channels 
denoted by c, o~ c is the sensitivity, C is the beam width, 
and l c = 2\/21n2/0 c is the corresponding "cut-off" mul- 
tipole. Eauation ll5l assumes that different channels pro- 
vide independent constraints. We follow previous the- 
oretical "error forecast" work in adopting this assump- 
tion; however, we note that this implicitly assumes that 
all foregrounds were perfectly removed from the tempera- 
ture and polarization maps. In reality, imperfect removal 
of foregrounds will induce correlations among the chan- 
nels, which will have to be taken into account in a refined 
analysis. Modeling foregrounds and the expected preci- 
sion with which they can be removed is beyond the scope 
of the present paper. 

In this paper we focus on the Planck survey for con- 
creteness. This survey will measure temperature and po- 
larization anisotropy in three frequency bands, namely 
100, 143 and 217 GHz, with fractional sky coverage of 
/sky ~ 0.8. The parameters for this experiment are 
listed in Table [I] (taken from [36J]). The various C/s 
are calculated up to £ max = 2, 000, as appropriate for 
Planck, using KINKFAST 37], a modified version of 
CMBFAST [38| tailored for time-varying w. 



C. Fiducial Cosmology 

The Fisher matrix formalism estimates how well a sur- 
vey can distinguish a fiducial model of the universe from 
other models. The results depend on the fiducial model 
itself. Throughout the paper, we take a 7-dimensional 
parameterization of a spatially-flat (Oj. = 0), cold dark 



D. Survey Parameters 

To determine the detection mass limit M m [ n (z) in 
Eq. 0, a mass-observable relation is needed. We con- 
sider three surveys for our analysis, a flux limited X-ray 
survey, such as DUO, an SZE survey that is similar to the 
one to be carried out with the South Pole Telescope, and 
a weak lensing survey that is similar to that planned for 
the LSST. For the X-ray and SZE surveys, we impose a 
minimum mass of 10 14 /i _1 Af Q (if M min (z) as computed 
below turns out to be less than lO 14 ft. -1 M0) since less 
massive halos correspond to small clusters or groups and 
are likely to depart from the scaling relations adopted 
here. We do not impose this lower bound for the weak 
lensing survey since, in principle, the dark matter halos 
of low-mass clusters or groups should still produce shear 
signals with a well-defined mass-shear relation. 

For the D UO-like X-ray survey, we adopt a bolometric 
flux-mass relation of the form: 



f x (z)4ndl = A x M^ E 2 (z)(l + z) 



0* 



(17) 



where f x in units of erg s _1 cm~ 2 is the bolometric flux 
limit, di, in units of Mpc is the luminosity distance, M200 
is the mass of the cluster, and H{z) = H E(z) is the 
Hubble _parameter at redshift z. Following Majumdar & 
Mohr |j, we adopt \og w (A x ) = -4.159, /3 X = 1.807, 
and 7 X = as fiducial values. We model the DUO ob- 
servations as a combination of a "wide" survey, covering 
a sky area of AD, = 6,000 deg 2 with a bolometric flux 
limit of f x > 1.75xl0 -14 erg s" 1 cm~ 2 (corresponding to 
f x > 7 x 10~ 14 erg s" 1 cm~ 2 in the 0.5 : 10 keV band); 



and a "deep" survey, spanning Aft =150 dcg 2 with a 
bolometric flux limit of f x > 2.25 x 1CP 14 erg s _1 cm -2 
(f x > 9 x 1CT 15 erg s _1 cm" 2 in the 0.5 : 10 keV 
band). With these parameters, for our fiducial cosmolog- 
ical model, the wide survey yields ~ 10, 000 clusters while 
the deep survey yields ~ 1,500 clusters. These numbers 
are consistent with the existing data on the log N — log S 
relation for clusters in soft X-ray bands [4(| and also 
with independent estimates for the total number of clus- 
ters expected to be detected by DUO [2lj. 

For an SPT-like SZE survey, we adopt an SZE flux- 
mass relation: 

f sz {z)d\ = f(v)hcuA sz M^EV\z){l + z)i", (18) 

where f sz in units of mJy is the observed flux decrement, 
(1a in units of Mpc is the angular diameter distance, f(f) 
is the frequency dependence of the SZE distortion, and 
/icm is the mass fraction of the intra-cluster medium. 
We model the SPT-like SZE survey as a flux-limited 
survey with f sz > 5 mJy at 150 GHz. While this is 
an oversimplification, the threshold value approximately 
represents the total flux decrement of the smallest cluster 
that can be detected at ha significance with SPT. [53 
We also adopt the fiducial parameters log 10 (A sz ) = 8.9, 
(3 SZ = 1.68, j sz = 0, and covering Aft = 4,000 deg 2 . We 
also assume /icm = 0.12. With these parameters, for our 
fiducial cosmological model, this survey yields ~ 20, 000 
clusters. 

Finally, for the LSST-like WL survey, we follow 
Hamana et al. |4l| to find a relation between the di- 
mensionless shear and halo mass, given by 



k g = a(9 G ) 



Af vir /(7rr 2 ) 



(19) 



Here kg is the average shear within a Gaussian filter of 
angular size 9g ; M V \ T and r s are the mass and scale radius 
of a cluster with an NFW density profile (r s = r vir /c n fw); 
zi is the redshift of the cluster; and d,A is the angular 
diameter distance to the cluster. The coefficient a is 
given by 



f™dx (x/x 2 G ) exp (-x 2 /x 2 G ) f(x) 

log(l + Cnfw) ~ C n fw/(1 + C n fw) 



(20) 



where x — 4>/9 s a dimensionless angular coordinate, and 
xg = 9g/9 s corresponds to the smoothing scale, with 
8 S = r s/dA(zi) denoting the angular scale radius. The di- 
mensionless surface density profile f(x) is given by equa- 
tion (7) in Hamana et al. |41|. Finally, the mean inverse 
critical surface mass density is given by 



5T/ = 



i 4ttG . . . J™ dz dn/dz (l- X (zi)/x(z)) 

-a(zi)x(zi) oo , 

J dz dn/dz 



where a is the scale factor and x denotes the comoving 
radial distance (valid for the flat universe with ftfc = 
we are assuming). 



a 

-a 

N 

-a 

-a 




FIG. 1: Expected number of clusters per redshift per unit 
solid angle for the fiducial cosmology. 



We model the LSST-like survey to have a constant 
detection threshold of kq = 4.5cr n oise |42|. The noise 
is given by the ratio of the mean ellipticity dispersion 
(cr e ) of galaxies and the total number of background 
galaxies contained within a smoothing aperture Qq |43| . 
°noisc = t7 e/(47T0G n s)- ^ e adopt <j e =0.3 and a number 
density of background galaxies n g — 65 arcmin -2 |44| . 
and the angular smoothing scale 6g = 1 arcmin, which 
corresponds to optimal S/N for the range of cluster 
masses and redshifts we considered. We take the survey 
to cover Ail = 18, 000 deg 2 and to extend over the clus- 
ter redshift range 0.1 < z < 1.4. This yields - 200,000 
clusters for our fiducial cosmology. 

Evidently, the various surveys considered here span dif- 
ferent redshift range, sky coverage and sensitivity. It is 
therefore useful to compare the expected number of clus- 
ters per redshift and unit solid angle for each survey for 
our fiducial cosmology. This is shown in Fig. Q] We see 
that the SPT-like survey is the most sensitive probe at 
high redshift, a consequence of the fact that the mass 
limit for SZE surveys is nearly redshift-independent. In 
comparison, the DUO X-ray and LSST-like surveys drop 
more sharply with redshift. Note that the fixed "mass 
floor" of \Q 1A h~ l Mq determines the number of clus- 
ters at z < 0.2,0.25,0.5 for DUO wide, SPT-like and 
DUO deep, respectively. In comparison, the LSST-like 
counts are dominated by small clusters or groups with 
M < 1O 14 /i- x M . 



III. RESULTS 

Before discussing our results for the cluster surveys, 
we first summarize constraints from the CMB alone. 
Our projections for the Planck satellite, using the seven- 
parameter Fisher matrix, are listed in Table [H] and are 
consistent with well-known previous forecasts [35|, EM 
I46J ]. This table shows the power of the CMB in con- 
straining the matter and baryon density, fl m h 2 and flbh 2 , 



TABLE II: Estimated Cosmological Parameter Constraints 
from Planck. 



Planck Survey 




A^DE 


0.035 


AQ. m h 2 


0.0012 


A(78 


0.041 


Aw 


0.32 


Aw„ 


1.0 


An b h 2 


0.00014 


An s 


0.0035 



TABLE III: Estimated Cosmological Parameter Constraints 
from DUO. The dN/dz column includes priors from WMAP: 
AQ. b h 2 = 0.0010, and An s = 0.040. 



Survey and Parameter Constraints 


dN/dz 


Pc(k) 


DUO Wide (6,000 deg 2 ) 






A^DE 


0.14 


0.037 


AQ m h 2 


0.25 


0.096 


Acrg 


0.16 


0.10 


Awo 


0.16 


0.59 


Aio„ 


0.92 


3.2 


AQ b h 2 


0.0010 


0.023 


An s 


0.040 


0.18 


DUO Deep (150 deg 2 ) 






A^DE 


0.097 


0.11 


Att m h 2 


0.33 


0.25 


Acrg 


0.040 


0.25 


Awo 


0.29 


0.78 


Aio„ 


2.5 


3.7 


AQ b h 2 


0.0010 


0.059 


An, 


0.040 


0.49 



respectively, as well as the spectral tilt, n s - However, as 
is well known, the equation of state of the dark energy, 
parameterized by wq and «v, is poorly constrained by 
CMB observations alone [471 |48| . This is because the 
dependence of the CMB power spectrum on these two 
parameters comes mainly from the distance to last scat- 
tering, c?LS) which involves a double integral of w(z): 



dhS 



dz 



3000 Mpc J ^/n m h2(l + z)3 + (l-n m )h2g(z) ' 
where z roc is the redshift of recombination, and where 



g(z) = exp\ 3 



(1 + w(z'))dz' 

1 + z' 



(23) 



Since wq and w a only appear inside this double integral, 
there is a severe degeneracy that can leave c^ls nearly 
invariant under changes in these parameters. 

Table IIIII summarizes the results for the DUO wide 
and deep surveys, both from cluster counts (dN/dz) and 



power spectrum (P c (k)) 58] . This table addresses the is- 
sue of the relative merits of survey size versus depth. 
Starting with the counts, we see that the constraints 
for the wide and the deep surveys are of the same or- 
der, even though the latter yields about 7 times fewer 
clusters. This is because the deep survey, despite its 
limited angular coverage, measures a higher fraction of 
high-redshift clusters. For the power spectrum, however, 
the constraints are most sensitive to the total number 
of clusters. Indeed, the erro rs on most p arameters dif- 
fer by roughly a factor of ^/A^idc/^dccp ~ 2.6, where 
TVwidc and ./Vdccp are the total number of clusters for the 
respective surveys. The power spectrum P c (k) delivers 
good constraints on the densities ($Ts) and on a% in the 
wide survey, but has little sensitivity to wq and w a . 

Table IIVI addresses the issue of the relative merits of 
dN/dz, P c (k), their combination, and their combination 
with CMB anisotropy data from Planck, for all three sur- 
veys. The top third of Table ITVI lists the results of com- 
bining DUO wide and deep. Of particular interest is the 
third column from the left, which shows the constraints 
obtained by adding the Fisher matrices for dN/dz and 
P c (k) for the combined survey. These are, in short, the 
most optimistic error bars from DUO alone. The ta- 
ble also illustrates the power of combining cluster counts 
with two-point function statistics. Indeed, the combined 
error bars (column 3) for J7de and fl m h 2 are about two 
times smaller than those derived from either dN/dz (col- 
umn 1) or P c (k) (column 2) alone. Finally, we see in the 
last column that combining DUO and Planck further re- 
duces the uncertainty on wq and w a by approximately a 
factor of two. This underscores the complementarity of 
cluster and CMB data in uncovering the nature of the 
dark energy. 

The middle third of Table IIVI shows our results for 
the SPT-like SZE survey. Overall, the constraints on 
the cosmological parameters are similar to those available 
from the DUO-like survey. 

The bottom third of Table IIVI shows our estimated 
parameter uncertainties for the LSST-like cluster sur- 
vey. Comparing the third column with the previous two 
again confirms the power of combining counts with power 
spectrum. The constraints on wq and w a are of the or- 
der of a few percent. These remarkably tight bounds 
(comparing favorably with those from the Planck survey 
[Table HT| for all cosmological parameters except fl m h 2 , 
tlbh 2 and n s ) are the result of the very high cluster yield 
of this survey. To examine the sensitivity of these results 
to the inclusion of the lowest mass clusters, we follow 
a more conservative approach by imposing a minimum 
mass of 2x 10 1 h~ l M@. This reduces the number of clus- 
ters to ~ 50, 000. The constraints on w(a) from LSST 
alone degrade by a factor of about two (consistent with 
y/~N scaling of statistical errors), to Awo = 0.030 and 
Aw a — 0.23. When combined with Planck, the errors 
are nearly unaffected: Awo = 0.024 and Au> a = 0.070. 
Therefore we are confident that with enough clusters and 
combining with Planck, LSST can constrain w a to a few 



TABLE IV: Estimated Cosmological Parameter Constraints from Clusters and CMB Combined. The dN/dz column includes 
priors from WMAP: AQ b h 2 = 0.0010, and An s = 0.040. 



Survey and Parameter Constraints 


dN/dz 


Pc(k) 


dN/dz + P c {k) 


dN/dz + P c (k) + Planck 


DUO Combined 










A^DE 


0.011 


0.032 


0.0074 


0.0064 


Afl ra ft 2 


0.022 


0.084 


0.0098 


0.00041 


A<7 8 


0.016 


0.088 


0.012 


0.011 


Aw 


0.10 


0.45 


0.096 


0.061 


Ato a 


0.48 


2.3 


0.45 


0.19 


AQ b h 2 


0.0010 


0.021 


0.0010 


0.00011 


An s 


0.040 


0.15 


0.033 


0.0024 


SPT-like Survey (4,000 deg 2 ) 










AfiDE 


0.036 


0.033 


0.014 


0.0097 


Attmh 2 


0.049 


0.056 


0.0083 


0.00027 


A<7 8 


0.031 


0.064 


0.018 


0.012 


Aw 


0.22 


0.41 


0.15 


0.082 


Aw a 


0.59 


1.8 


0.46 


0.18 


An b h 2 


0.0010 


0.014 


0.00099 


0.00011 


An s 


0.040 


0.094 


0.029 


0.0023 


LSST-like Survey (18,000 deg 2 ) 










A^DE 


0.0053 


0.0080 


0.0024 


0.0023 


Afl ra ft 2 


0.026 


0.021 


0.0048 


0.00024 


Aag 


0.0035 


0.022 


0.0025 


0.0024 


Awo 


0.051 


0.10 


0.024 


0.023 


Aw a 


0.086 


0.47 


0.077 


0.061 


AQ. b h 2 


0.0010 


0.0050 


0.00097 


0.00010 


An s 


0.040 


0.040 


0.015 


0.0022 



percent. Finally, we find that it is essential to include the 
cosmology-dependence of the limiting mass for the LSST 
survey. Repeating our analysis adopting the redshift- 
dependent mass limit from the fiducial cosmology, and 
not allowing it to vary with cosmology, results in an in- 
crease by a factor of 3-4 in the uncertainties. 

In Table E] we repeat the analysis for the DUO and 
SPT-like surveys, but this time taking into account the 
uncertainty in the structure and evolution of clusters. In 
other words, we require that the cluster surveys constrain 
not only the cosmology, but also the parameters of the 
mass-observable relation. For DUO, this is modeled by 
including the parameters A x , (3 X and ^ x of Eq. (fTTjl in the 
Fisher matrix analysis. Similarly, for the SPT-like sur- 
vey, we include A sz , f3 sz and ^ sz from Eq. (fTHfl . For both 
surveys, we see that adding self-calibration increases the 
error on wq and w a by a factor of < 2 in comparison 
with the corresponding results in Table llVl Overall, we 
see that including self-calibration parameters still yields 
very good constraints on the cosmology. The constraints 
from dN/dz+P c (k) and dN/dz+P c (k)+Pla,nck are shown 
graphically in Figs. El and including self-calibration for 
the DUO-like and SPT-like surveys. 

Coming back to Table IIVI and comparing the first col- 
umn from the left (dN/dz) with the third (dN/dz + 



P c (k)), we note that, for DUO, adding the power spec- 
trum does not significantly tighten the constraints on wq 
and w a in this case. However, for the "self-calibration" 
case (Table 0, combining the two methods helps greatly 
on all the constraints. For SPT-like and LSST-like sur- 
veys, combining two methods always gives stronger con- 
straints. This is clearly illustrated in Fig.|3 

An obvious method to cross-check for systematic ef- 
fects due to cluster structure and evolution, and to im- 
prove constraints, is to combine the X-ray and SZE data. 
We have found that the uncertainties in wq and w a reduce 
to 0.083 and 0.35 when the self-calibrated SPT and DUO 
samples are considered in combination, and to 0.060 and 
0.22 when Planck is added to this combined sample. 

Uncertainties in cluster structure and evolution should 
be less severe for WL signatures, which probe the dark 
matter potential directly. However, to compare the WL 
error forecasts more fairly with the SZ and X-ray pre- 
dictions, in Table IVTl we show LSST predictions that in- 
corporate two additional uncertainties inherent to WL 
selection: false detections and incompleteness. We relied 
on published results of numerical simulations |4lt |42J | to 
calibrate our results: we assumed 25% of detections are 
false, and 30% of real clusters are undetected. These ef- 
fects increase the parameter errors, relative to those in 



TABLE V: Parameter Constraints including Self-Calibration. The dN/dz column includes priors from WMAP: AQ, b h 2 
0.0010, and An s = 0.040. 



Survey and Parameter Constraints 


dN/dz 


Pc{k) 


dN/dz + P c (k) 


dN/dz + P c (k) + Planck 


DUO Combined 










AQde 


0.030 


0.043 


0.015 


0.012 


Afl m h 2 


0.14 


0.091 


0.0098 


0.00067 


Acrg 


0.058 


0.12 


0.016 


0.013 


Awo 


0.44 


0.53 


0.20 


0.15 


A* 


1.2 


2.5 


0.82 


0.46 


AQ b h 2 


0.0010 


0.022 


0.0010 


0.00011 


An s 


0.040 


0.18 


0.034 


0.0027 


A log A x 


0.27 


0.18 


0.050 


0.037 


Ab x 


0.29 


0.18 


0.030 


0.027 


Ay x 


0.67 


0.63 


0.17 


0.081 


SPT-like survey 










Afi DE 


0.13 


0.037 


0.017 


0.014 


An m h 2 


0.65 


0.077 


0.0086 


0.00063 


Acrg 


0.14 


0.099 


0.019 


0.017 


Awo 


0.42 


0.46 


0.19 


0.14 


Aw a 


2.4 


1.8 


0.81 


0.40 


AQ b h 2 


0.0010 


0.018 


0.0010 


0.00011 


An s 


0.040 


0.15 


0.032 


0.0026 


A log A sz 


0.59 


0.35 


0.12 


0.056 


Ab sz 


0.79 


0.35 


0.12 


0.063 


A^ sz 


1.7 


0.62 


0.19 


0.057 



TABLE VI: Calibrated Cosmological Parameter Constraints from LSST and CMB Combined. The dN/dz column includes 
priors from WMAP: AQ, b h 2 = 0.0010, and An s = 0.040. 



Parameter Constraints 


dN/dz 


Poik) 


dN/dz + P c {k) 


dN/dz + P a (k) + Planck 


LSST-like Survey 










AfiDE 


0.0081 


0.012 


0.0037 


0.0033 


An m h 2 


0.038 


0.032 


0.0059 


0.00024 


Acrg 


0.0054 


0.034 


0.0038 


0.0037 


Aw 


0.079 


0.16 


0.037 


0.036 


AlOa 


0.13 


0.72 


0.12 


0.093 


AQ b h 2 


0.0010 


0.0077 


0.00098 


0.00010 


An s 


0.040 


0.062 


0.021 


0.0022 



the bottom third of Table [El by « 50% (except for fl b h 2 
and n s , which are significantly determined by the WMAP 
priors and are less affected). See the next section for a 
more detailed discussion. 



IV. DISCUSSION 

In the previous section, we derived constraints on cos- 
mological parameters from future SZE, X-ray, and WL 
surveys. Our Fisher matrix approach should be inter- 
preted as yielding lower limits on the achievable statisti- 
cal errors. In our analysis, we adopted unique relations 
between the observables and cluster mass, which come 



from simple models of cluster structure and evolution. 
We required the SZ and X-ray surveys to "self-calibrate" 
and constrain structure and evolution parameters simul- 
taneously with cosmology. However, the power-law form 
of the relations we adopted will have to be further con- 
strained. This should be feasible by combining the three 
observables in the different wave bands for a subset of 
the samples, and by adding new observables (such as the 
shape of the mass function, the angular size, velocity dis- 
persion) which we have not considered here. In the case 
of the WL sample, we relied on results from numerical 
simulations to calibrate the mass-observable relation, an 
approach that can be refined with a larger suite of simu- 
lations in the future. 
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FIG. 2: Constraints on dark energy parameters: wo and w a (left), wo and JIde (right) for a DUO-like X-ray survey (top), 
an SPT-like SZE survey (middle), and an LSST-like weak lensing survey (bottom). The three curves in each figure show the 
constraints available from dN/dz (dashed), P c (k) (dotted), and from their combination (solid). The star-shaped symbol at the 
center of each figure indicates our fiducial cosmology. The constraints for X-ray survey and SZE survey are calculated for the 
self-calibration case. The constraints are marginalized over all other cosmological and relevant structure parameters. In all 
cases, the constraints from the combination of dN/dz and P c (k) are at least a factor of two stronger than from either method 
alone. 



The Fisher matrix technique, especially the way of 
combining independent constraints by summing individ- 
ual Fisher matrices, allows us to explore the physical ori- 
gin of the cosmological information. A glance at Eqs. JSJ 
and IjlUI) reveals that cosmology enters through several 
physical quantities into the Fisher matrices, such as the 
cosmic volume, growth function, transfer function, bias, 
etc. Unfortunately, the marginalization over all other 



parameters involves a nonlinear inversion of the Fisher 
matrix, which makes isolating the various sources of in- 
formation difficult. As an example, we repeated the anal- 
ysis above, but keeping either the volume factor, fc^/c||, 
or the bias parameter at their values in the fiducial model 
(i.e., excluding their derivatives from computing power 
spectrum Fisher matrices). We find that as a result, the 
constraint on some of the parameters improve, while oth- 
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FIG. 3: Constraints on dark energy parameters: wq and w a (left), wo and JIde (right), by combining a DUO-like X-ray survey 
(top), an SPT-like SZE survey (middle), or an LSST-like weak lensing survey (bottom) with Planck-like CMB observations. 
The three curves in each figure show the constraints available from the combination of dN/dz and P c (k) (dotted), CMB alone 
(dot-dashed), and from the combination of all of the above (solid). As in Fig. |^1 the constraints for X-ray survey and SZE 
surveys are calculated for the self-calibration case, and constraints are marginalized over all other relevant parameters. Note 
the different scales of the horizontal and vertical axes compared to Fig. |5] 



ers degrade, which does not offer a useful description of 
the amount of information the volume factor or the bias 
parameter provides. 

In addition, we explored the implicit assumption made 
above that the bias parameter b can be precisely modeled. 
We take a similar approach to that used in addressing 
"self-calibration", by modeling the bias as b = bo(z)(l + 
z) lb . We effectively include an additional new parameter, 
7b, in the Fisher matrix analysis. Note that a constant 



factor of normalization would be degenerate with as ■ In 
the fiducial model, bo(z) follows from Eq. (J5J, and 7b = 0. 
The constraints from cluster power spectrum including 
the "non-standard evolution parameter" jf, causes only 
a minor degradation (under 10%) of the constraints on w 
and w a (relative to column 2 in Table llVJI . suggesting that 
the bias factor did not drive the cosmological constraints 
we derived from the power spectrum. 

The results from the power spectrum also depend on 
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the assumption of the extent of the linear regime. To 
quantify the importance of the small-scale modes, for 
the example of the SPT-like survey, we decreased £; max 
from 0.15 to 0.10 Mpc -1 . We found that this degrades 
the error bars by a factor of up to *=* 1.4. Similarly, 
decreasing fc max further to 0.075 Mpc^ 1 degrades all the 
error bars by a factor of up to w 2. (See also 14] for a 
detailed discussion of the dependence of the constraints 
separately on fcx.max and fc|| iinax .) 

In Table IIVI we have assumed that the mass- 
observable relations apply exactly. However, for real sur- 
veys, the mass-observable relation is likely to have a non- 
negligible scatter. We estimated the magnitude of scatter 
that would bias the inferred parameters by an amount 
comparable to their statistical errors, for the example of 
the SPT survey [27|. We assumed that at a fixed flux 
fsz, the mass Miso has a Gaussian distribution, with a 
mean given by Ea.ll8lfbut converted from M200 to Miso), 
and a fractional r.m.s. deviation of SM/M — <tm- We 
then re-computed the cluster abundance dN/dz (which 
is increased in the presence of scatter). We found that 
at <7j\/ ss 7%, the change caused in the total number of 
detectable clusters is comparable to the y/~N error that 
is used to obtain the constraints on cosmological param- 
eters. A scatter larger than ~ 7% would therefore be 
important, and it would have to be modeled (i.e. by 
parameterizing the scatter, calibrated directly from ob- 
servations) in the analysis of real data. 

A further complication, likely most relevant for the WL 
survey, is the presence of false detections. In our analysis 
of the LSST-like survey, we used a constant shear S/N 
to select clusters, which directly measures how well the 
shear signal stands out above the ellipticity shot noise 
of the background galaxies. False detections can result 
from statistical fluctuations in these ellipticities. Fur- 
thermore, physical structures, projected along the line of 
sight to a given cluster, constitute additional background 
noise, and can result in a false detection of two or more 
mass concentrations (projected along the line of sight but 
physically separated in redshift) . Hennawi & Spergel |43 
have done a comprehensive study of mass-selected clus- 
ters using A-body simulations (see also |4lll49J for other 
numerical studies of WL cluster selection). Defining ef- 
ficiency as the fraction of the peaks in the mass map 
that corresponds to real clusters, and completeness as 
the fraction of real clusters we can detect, they find in 
their simulation that for a 4.5 standard deviation detec- 
tion threshold, sa 75% of the detected clusters are real 
and ps 70% of the clusters can be detected |41|. They 
emphasize that the expected cluster distribution, includ- 
ing false detections and missing clusters, can be reliably 
calculated for any cosmological model since the simu- 
lations depend only on gravity. Thus, false detections 
can be subtracted, and their presence serves only to in- 
crease the statistical error. These effects would increase 
the parameter errors in the bottom third of Table llVl by 
~ 50% (summing in quadrature the errors on the num- 
ber of real clusters and on the number of false-detections, 



^/((1/e - 1) + l/e)/c - 1 sw 50%, where e = 75% is the 
efficiency and c = 70% is the completeness as we defined 
above). 

We note that the study in 42] utilized only the "tomo- 
graphic redshifts" derived from the shear map itself (us- 
ing photometric redshifts only for the background galax- 
ies), which have large uncertainties. However, it should 
be feasible to routinely measure photometric redshifts of 
the cluster member galaxies, as well, and this additional 
information could be used to reduce the rate of false de- 
tections. 

Obtaining redshifts is a general issue relevant for all 
three future surveys. While accurate (spectroscopic) red- 
shifts are not required for either the dN/dz or P(k) tests, 
photometric redshift estimates will be required to utilize 
the distribution of clusters in z, and also to use high-fc|| 
modes of the power spectrum. In the case of the LSST 
survey, we adopted a maximum redshift z max = 1-4, the 
redshift out to which it should be feasible to obtain pho- 
tometric redshifts with the planned filters and sensitiv- 
ities I23. In the case of DUO, most of the clusters in 



the wide survey ( 



> 



will have photometric redshifts 



from the overlapping Sloan Digital Sky Survey (SDSS) 
galaxies; the redshifts of the more distant galaxies in the 
deep survey will be obtained in optical follow-up pro- 
grams. Obtaining photometric redshifts of a large frac- 
tion of clusters in the SPT survey beyond z > 1.4 will be 
challenging. However, we have re-computed our results 
for SPT ignoring all clusters beyond z = 1.4, and found 
that none of the constraints degraded by more than a few 
percent. 

It has often been argued that weak lensing only in- 
volves gravitational deflection of light, so the mass-shear 
relation can be determined without fully understanding 
the complex baryonic physics. Clearly this is an ideal- 
ization, since baryons can cool and contract to establish 
a more centrally-condensed density profile than the dark 
matter, altering the total gravitational field. To quan- 
tify the resulting change in the WL shear, we followed a 
simple spherical model in [50j, in which a fixed fraction 
/cool of baryons inside a 'cooling radius' are allowed to 
cool and condense by a further factor of ~ 10 after viri- 
alization. The effect of such a "cooling flow" on the mass 
profiles of halos is relatively small, except in the inner 
regions. In models with / coo ] up to 30 %, we find that 
the change in the Gaussian-averaged shear kq can be up 
to 15%. This reflects the enhanced effect of baryons in 
the core of the cluster, which is emphasized by the Gaus- 
sian weighting. However, this bias only affects the initial 
selection of clusters. In practice, one will go back and 
extract the total mass of the cluster from the shear map 
in a separate analysis [51|. We find that the total mass 
enclosed within the virial radius is changed by a much 
smaller amount, < 1%. 

In this work, the constraints on cosmological parame- 
ters were derived assuming a fiducial model with a cos- 
mological constant as the dark energy, i.e., wq — — 1 and 
w a = 0. As mentioned above, the choice of the fiducial 
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model can, in principle, strongly affect the predicted er- 
rors. To assess the robustness of our derived errors, we 
repeated the analysis for a fiducial model with a strongly 
time-dependent equation of state, namely wq = —0.8 and 
w a = 0.3, which is also consistent with current data. We 
find in this case that all error bars derived from dN/dz 
and P(k) remain essentially unchanged, except for Aw a , 
which shrinks by about 25%. Therefore, by choosing a 
pure A model as our fiducial cosmology, we obtain a con- 
servative estimate of the error on w„ . 



V. CONCLUSIONS 

Three main lessons can be drawn from our work. First, 
the above results provide further support for the useful- 
ness of clusters as an independent and important probe 
of the dark energy. At a general level, the main advan- 
tage of clusters over other measurements lies in the fact 
that the physics involved is predominantly gravitational 
and hence simple; moreover, the observations are at low 
redshift and therefore straddle the epoch of cosmic accel- 
eration. 

We have seen that even cluster counts (dN/dz) alone 
can yield comparable or even better constraints on wq 
and w a than a full-sky, high-precision CMB survey such 
as Planck. Our constraints our summarized in Table llVl 
Of course, this assumes perfect knowledge of the mass- 
observable relation, i.e., that the clusters are standard 
candles. We have shown, however, that even after al- 
lowing for uncertainty in cluster evolution and structure, 
and using the survey itself as a calibration tool, the er- 
rors on Wq and w a still compare respectably well with 
Planck, as seen from the third column of Table ]V\ The 
constraints from cluster surveys also compare favorably 
with other cosmological probes. For instance, current 
constraints from type la supernovae are on the order of 
Awo ta 0.3 and Aw a ss 1.6 J52J, assuming a strong prior 
of Afi m = 0.04 from other measurements. The expected 
errors from the Supernovae Acceleration Probe (SNAP) 
are Awq_~ 0.08 and Aw a ~ 0.3, with the same prior 
on fl m |53j |. Relaxing this prior results in errors on the 
order of Awq ss 0.2 and Aw a Ril.O [47|, comparable to 
our projected constraints from DUO or SPT. 

A second lesson we wish to draw is the power of com- 



bining different cluster observablcs in constraining cos- 
mological parameters, in this case cluster counts (dN/dz) 
and two-point statistics (P c (k)). The most striking illus- 
trations of this are the LSST-like survey (see Table IIV|) 
where Awq is a factor of ~ 2 smaller after combining 
the two methods; and the SPT-like self-calibrated case 
(see Table where both Awo and Aw a shrink by a fac- 
tor of > 2. While we have focused here on dN/dz and 
P c (k), other cluster observables that could be included 
are three-point correlation functions, assuming a large 
enough survey, information from the shape of the mass 
function, dN/dM, as described in [ll|, or scaling rela- 
tions [54( . We leave a more comprehensive study of such 
additional cluster observables for future work. 

Finally, our work underscores the complementarity of 
cluster observables to CMB anisotropy. In general, we 
have seen that adding the information from a Planck- 
like survey reduces Awq and Aw a by a factor of 2 or 
so (although the LSST-like WL survey has already very 
small errors, and adding Planck data improves the con- 
straints by a more modest amount). Thus clusters not 
only constitute a powerful cosmological probe in their 
own right, but also help in alleviating CMB observations 
from some of their well-known degeneracies. 

Future work includes the addition of other cluster ob- 
servables, as mentioned earlier, as well as using clus- 
ter surveys to place useful constraints on the neutrino 
mass 1551. 
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